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Abstract 

We report robust calculations of various low energy QCD hadronic prop- 
erties. We use a multi-rank separable expansion for the gluon propagator 
which greatly facilitates the numerical computations within the Global 
Colour Model for QCD. The parameters for the propagators are deter- 
mined by fitting experimental values for f n and the ir and a± meson 
masses. 

PACS numbers: 12.38.Lg, 13.75. Cs, ll.10.St, 12.38.Aw 



The computation of the low energy properties of QCD is a difficult non-perturbative 
problem in quantum field theory. Here we report new results using the Global Colour 
Model (GCM) jjj approximation to QCD. These robust computations are made feasible 
by the use of multi-rank separable expansions of the gluon propagator. While good 
progress |2| has been made in computing the gluon propagator from first principles in 
QCD, here we adopt the procedure of determining an effective gluon propagator by 
using a convenient separable form, and selecting the parameters therein by fitting f n 
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and the n and a± meson masses to experimental values. The validity of both the GCM 
and the effective separable-form gluon propagator are demonstrated by then computing 
numerous other low energy hadronic properties. A non-separable translation-invariant 
form for the effective gluon propagator, appropriate to low energy hadronic states, is 
then reconstructed from the separable form. 

An overview and an insight into the nature of the non-perturbative low energy 
hadronic regime of QCD is provided by the functional integral hadronization of QCD 
||, 0]. This amounts to a dynamically determined change of functional integration 
variables, from quarks and gluons, to bare hadrons 

/ VqT>qVAexp(-S qcd [A, q, q]+rjq + qrf) « 

j VTTVNVN...exp(-S had [K,...,N,N,..] + J 7r [r},ri]TT + ..) (1) 

The final functional integration over the hadrons gives the hadronic observables, 
and amounts to dressing each hadron by, mainly, lighter mesons. This functional in- 
tegral transformation cannot yet be done exactly. The basic insight is that the quark- 
gluon dynamics, on the LHS of (1), is fluctuation dominated, whereas the RHS is not, 
and for example the meson dressing of bare hadrons is known to be almost pertur- 
bative. In performing the change of variables essentially normal mode techniques are 
used H . In practice this requires detailed numerical computation of the gluon propa- 
gator, quark propagators, and meson and baryon propagators. The mass-shell states 
of the latter are determined by covariant Bethe-Salpeter and Faddeev equations. The 
Faddeev computations are made feasible by using the diquark correlation propagators, 
which must also be determined. 

The first and easiest formal transformation results from doing the gluon integra- 
tions, leaving an action for quarks of the form 

S[q,q] = J q{ x ){- 1 . d +M)q{x)+ l - J j -( x )f u (y)D^x-y) + ^ J f^D^ ... (2) 

where j^(x) = q(x)^- r y fM q(x). The GCM is a model field theory for QCD based on a 
truncation of S\g, q] in which the higher order n-point (n > 3) functions are neglected, 
and only the gluon 2-point function D fll ,(x — y) is retained. 

The GCM is thus a quantum field theory that can also be considered to be defined 
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by the action 

[q(x)(--f.d + M + iA«—^)6(x - y)q(y)+ 

+ l -Al{x)D^{ l d)5{x-y)Al{y)^ (3) 

where the matrix (p) is the inverse of (p) , which in turn is the Fourier transform 
of Df JjU (x). This action has a global colour symmetry. The GCM is thus analogous 
to QED except for colour currents and the non-quadratic phenomenological form for 
D~„(p) in the pure gluon sector. The main purpose of this work is to report a robust 
and effective procedure for determining D^ v {p) and to demonstrate its general validity 
for a range of hadronic phenomena. This procedure is to use the separable expansion 
technique which proved very effective in the 1960's in studying non-relativistic few 
particle systems. 

Having made the GCM truncation in (2) it is possible to proceed further and to 
transform Q the quark functional integrations into the hadronic functional integrations, 
as in (1). If the additional approximation D^ u (x—y) — > g5 flu S(x—y) is made in (2), i.e. a 
contact coupling of the quark currents, then the NJL type models are obtained. If in (1) 
a derivative expansion of the complete non-local hadronic effective action is performed, 
then the Chiral Perturbation Theory (CPT) phenomenology is obtained. However 
in the GCM, with appropriate D^ u (x), all computations are finite and no cutoffs or 
renormalisation procedures are used. As well, using a mean field approximation, the 
soliton phenomenology for the baryons may be derived Q, and has been studied in ||. 

As shown here and elsewhere this approximation is surprisingly effective in de- 
scribing the low energy QCD determined hadronic properties, and suggests that some 
particular mechanism is responsible for its success. There are indications that the ne- 
glected terms in (2) may only be significant in those bound states which are not colour 
singlets || . In colour singlet states some colour neutrality may render the higher order 
terms ineffective. In the GCM the remaining gluon 2-point function is now regarded 
as an effective 2-point function: D fiu (x — y) — > D fll/ (x — y) e ff- There are then two 
key steps in the GCM: (1) the determination of this effective 2-point function, and the 
demonstration that the same function does well in determining a variety of hadronic 
observables, so that in some sense it is universal, and (2) the comparison of -D e // with 
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the true one obtained from QCD. This would establish to what extent the effects of the 
higher order terms have been accounted for by -D e //- From now on we shall call this 
D. 

Having specified D(p) the usual procedure would be to first determine the (con- 
stituent) quark propagators, in which m is the quark current mass, 

G{q) = (iA(q)q.j + B{q) + m)" 1 = -iq.-fa v (q) + a s (q) (4) 

using the non-linear Dyson- Schwinger equations (DSE), in Euclidean metric, 
T3( \ 16 f A n , v B(q)+m 



(2tt) 4 u HJ 'q 2 A(q) 2 + {B{q)+m) 2 ' 

M>> " 111,2 " I / Wf q * DI * ~ + + (6) 

where (5), using the convolution theorem for Fourier transforms, has a particularly 
simple form in coordinate space 

B(x) = ^D(x)a s (x). (7) 

Eqn.(5) is the more sensitive one of (5) and (6): the meson form factors are closely 
related to the extent of B(p), whereas A(p) is changing more slowly. For simplicity 
we have used a Feynman-like gauge, and the perturbative quark-gluon vertex function. 
One can choose to use the GCM in other gauges. If one chooses a Landau gauge form 

d,Ap) = (V - ^) a (p 2 ) 

then from the 0(4) invariance of B(p), we obtain that D(p 2 ) = |A(p 2 ) in (5). A similar 
effect occurs in the BSE because for the low mass states, with the confining quark 
propagator, the form factors are almost 0(4) invariant wrt the relative momentum 
dependence ||. 

The DSE only have the forms in (5) and (6) for a translation invariant action, as 
in (2) or (3). To solve these equations for various D(p) and to analytically continue 
A(s) and B(s) into the complex s-plane, where s = q 2 , when solving for meson, diquark 
and baryon propagators is particularly difficult. Hence we have studied the well known 
separable expansion technique. This greatly facilitates the solutions of the SDE, the 
BSE and the Faddeev equations. 



We first imagine expanding D(p — q) into 0(4) hyperspherical harmonics 

D{p -q) = D (p 2 , q 2 ) + q.pD^p 2 , q 2 ) + ... (8) 

where, for example, 

D (p 2 ,q 2 ) = - r df3 sin 2 f3 D(p 2 + q 2 - 2pqcosf3). (9) 

Jo 

However note that putting q = in (9) gives D{p) = Dq(p 2 , 0). Thus the full D(p) may 
be easily reconstructed from the first term in the RHS of (8). Note also that only the 
first two terms in (8) are needed in (5) and (6); higher order terms do not contribute 
as they are orthogonal to the measures therein. 

To facilitate the numerical computations we then introduce multi-rank separable 
expansions for each term 

D (p 2 ,q 2 )= E r ^ 2 ) r *(9 2 )- (10) 

i=l,n 

Introduction of the separable expansion clearly breaks translational invariance and 
must be regarded purely as a numerical procedure, much like a lattice breaks translation 
invariance. Translation invariance is restored as the rank of the separability is increased. 
The infrared hadronic region appears to be well described by a rank n = 2 form for Dq. 
The DSE then have solutions of the form 

B(s)=Y,b i r i (s),... (11) 

Equations for the 6j, .. are easily determined by substituting these forms back into (5) 
and (6), giving coupled transcendental non- linear equations easily solved by iteration. 
Then a s and a v are seen to have the form of sums 

a s( s ) = X! a s{s)i, a v (s) = a v(s)i, (12) 

i=l,n i=l, fe 

However this, in principle, procedure suffers from the defect that since D is un- 
known then some ad hoc choice of the parametrised Tj must be made. The resulting 
cr s (s)i and a v {s)i develop spurious singularities which impede the use of the quark 
propagator in meson and baryon computations. 

A much more robust and physically sensible procedure is to specify a parametrised 
form for the a s (s)i and a v (s)i of the chiral- limit quark propagator with known ana- 
lyticity properties. We then use the DSE in an inverse manner to compute the Tj(s). 
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By using entire functions the much speculated but as yet unproven quark confinement 
property can be ensured. Suitable forms are based on simple model solutions of the 
DSE 0]. These forms are 

, , f , x f , 2s-f3 2 (l-exp(-2s/(3 2 )) 
o- s {s)i = Ciexp{-diS), a v {s) = , (13) 

where a v (s) has only a rank k = 1 expansion here. 

From these parametrised forms we can easily determine the various parameters 
{bi,..}, which will depend on the basic, m = 0, chiral-limit parameter set Vq = 
{ci, C2, d%, d,2, f3}- For example from (5) we obtain 

* = T*F** s KaJ%w (14> 

in which B(s) = B(s)i + B(s)2, and 

B(s)i = a s (s)i/(scj v (s) 2 + a s (s) 2 ). (15) 

The basic procedure is to find those values of these parameters Vq which allow the 
fitting of some small set of hadronic observables to the experimental values. We have 
chosen f n (see [pfl for the expression), which probes exclusively the space- like region 
of the quark propagators, and the a\ meson mass which extends the probe into the 
time-like region. The a\ meson is chosen because all two-quark states considered here, 
except for the 1~ diquark correlation, have lower mass. Thus the region in the complex 
s-plane where the quark propagator is needed has been probed by our fitting procedure. 
A space-like only fitting procedure would require unreliable and untested extrapolations 
into the time-like region when computing other hadronic observables. The tt mass is 
needed to mainly determine the averaged u and d current masses. 

The a± and tt masses are computed using the Bethe-Salpeter equations (BSE), 
keeping only the dominant Lorentz amplitude ||. The BSE are solved with the loop 
integration in the Euclidean metric, while the meson momentum, in the rest frame 
P = (0,iM), is in the time-like region of the Minkowski metric. This mixed metric 
is typical of hadron calculations in the GCM, and ensures that the quark propagators 
stay as close as possible to the real positive s-axis, which is the region most studied for 
the quark and gluon propagators. As is well known a separable expansion for the kernel 
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of an integral equation, here the gluon propagator in the BSE, reduces that equation 
to algebraic form. The solution of the BSE is thus rendered almost trivial. 

We now consider the inclusion of quark current masses in the computations. For 
any set of Vo values, which implicitly parametrise the gluon propagator, we can compute 
the non-chiral quark propagator by returning to solve (5) and (6) with m ^ 0. The 
resulting a s and a v now depend on m. However the approximate DSE is only robust 
in the space-like region. Hence we have fitted, in this region, a s {s;m) and a v {s;m) to 
the forms in (13). This ensures that the confinement ansatz continues to hold. It is 
this non-chiral quark propagator which is needed in the 7r and a\ mass computations. 
The value of the average u/d quark current mass m is varied, along with the Vo values, 
during the fitting procedure. Once the Vo values were known we were able to use the 
DSE for the strange quark. Various values of the strange quark current mass were used 
until the K-meson BSE produced the experimental K meson mass. 

The determined parameter values Vo are shown in Table 1, and the values of the 
fitted observables are shown in Table 2. These three observables are sufficient to give a 
robust determination of the parameter set and the u/d quark mass. Table 2 also shows 
some of the predicted observables which test the universality of the quark and gluon 
propagators. Figs. 1 and 2 show the form of the chiral-limit a s and a v . 

Now we consider the task of reconstructing the translation invariant form for the 
effective gluon propagator. We can easily do this using D(p) = Do(p 2 ,0) or (7) for 
D(x). Hence if the translation invariant constituent quark propagator is known, and 
the separability technique does not compromise that property, then the full D(p) or its 
Fourier transform D(x), determined finally by the values of the parameter set Vo, may 
be constructed. It is this D(p), but with the GCM analysis done in the appropriate 
gauge, which should be compared with future direct computations of the gluon propa- 
gator. Using the separable technique an explicit expression for D(p) may be obtained. 
We have 

D{p) = D (p 2 , 0)=J2 r^ 2 )F;(0) = ^B(p%B(0) t 

% i * 

_ v 1 a s (0)i a s {p 2 )i 
an expression which becomes increasingly more accurate as the rank of the separable 
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expansion is increased. With the parameter set in Table 1, (14) gives b\ = 0.02672 GeV 2 
and 62 = 0.02395 GeV 2 , and the resulting D(p) is shown in Fig. 3. Note that the meson 
data mandates that the gluon propagator has a strong IR component (corresponding 
to a large distance confining effect) and a longer range component (corresponding to a 
medium distance effect), as also seen in the QCD studies in ||. 

The predicted hadronic observables in Table 2 are mainly self-evident. The details 
of the GCM computational techniques have all been discussed in the literature. This 
set of hadronic predictions is obtained without cutoffs or renormalisation procedures. 
The meson and diquark masses are from BSE computations, while the nucleon-core 
mass (equivalent to the quenched approximation in lattice QCD) is from a covariant 
Faddeev computation H 10 1 keeping only the + diquark correlation. The separable 
expansion for the gluon propagator leads to a separable form for the diquark correla- 
tion propagator. The original nucleon-core Faddeev computations (!(]] assumed such a 
separable form, see || for details. 



Later computations of the nucleon core frTl| , 12, 13, 14] usually had free parameters 
that are adjusted to give the experimental nucleon mass. A full computation of the 
nucleon mass in the GCM requires the inclusion of the 1 + diquark correlation, and the 
dressing of this nucleon core by mesons (equivalent to including quark loops in lattice 
models). 

The constituent quark masses arise from the most probable value of the quark 
running masses (B(s,m) + m)/A(s) in the BSE kernel integrations, which is at s ~ 
0.3GeV 2 . Expressions for the tt — ir scattering lengths and the pion charge radius r w 
(without re-scattering corrections) are given in In [|J the MIT bag phenomenology 
was derived from the GCM using a mean field soliton approximation. This gave an 
expression for the MIT bag constant in terms of B(s) and A(s). 

The good values of the predicted observables suggest that the GCM is capable of 
accurately describing low energy QCD, and that the gluon propagator multi-rank sep- 
arable expansion is a particularly useful and robust computational tool. By including 
more experimental data in the fitted observables larger rank expansions could be used. 
This would result in a more accurate mapping out of the gluon and constituent quark 
propagators in the momentum range appropriate to low energy hadronic physics. 
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RTC acknowledges useful conversations with C.J. Burden, CD. Roberts, P.C. 
Tandy and M. Thomson. 

Addendum: Frank and Roberts (T^j have reported model gluon propagator results 
also from a fit to meson data using the separable expansion technique. Burden et al. 
|20t| have also used the separable expansion technique for the gluon propagator, but in 
particular they have studied the meson BSE equation beyond the dominant Lorentz 
amplitude approximation. The results of both works should be compared with the 
results reported here. 
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Figure 1: Plot of extracted a s (p) {GeV~ l ) plotted against p 2 (GeV 2 ). 

Figure 2: Plot of extracted a v (p) plotted against p 2 (GeV 2 ). 

Figure 3: Plot of extracted D(p) (GeV~ 2 ) plotted against p 2 (GeV 2 ). 



Table 1: Quark Propagator Parameters: Vq 

ci 0.5200G6V- 1 c 2 1.1794GCV- 1 

di 2.0737GeV- 2 d 2 4.7214GeV" 2 

(3 0.5082GeV 
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Table 2: Hadronic Observables 



Observable 



Theory Expt. /Theory [Ref] 



Fitted observables 

u 

a\ meson mass 
7r meson mass 

K meson mass (for m s only) 



93.00MeV 
1230MeV 
138.5MeV 
496MeV 



93.00MeV 
1230MeV 
138.5MeV 
496MeV 



Predicted observables 

(m u + m d )/2 
m s 

u> meson mass 
a(] 7r — 7r scattering length 
Oq 7r — 7r scattering length 
a\ 7r — 7r scattering length 
a® 7r — 7r scattering length 
a| 7r — 7r scattering length 
TV pion charge radius 
nucleon-core mass 
constituent quark rms size 
chiral quark constituent mass 
u/d quark constituent mass 
s quark constituent mass 
+ diquark rms size 
+ diquark constituent mass 
1 + diquark constituent mass 
0~ diquark constituent mass 
1~ diquark constituent mass 
MIT bag constant 

MIT nucleon-core mass (no cm corr. 



6.5MeV 
135MeV 
804MeV 
0.1634 
-0.0466 
0.0358 
0.0017 
-0.0005 
0.55fm 
1390MeV 
0.59fm 
270MeV 
300MeV 
525MeV 
0.78fm 
692MeV 
1022MeV 
1079MeV 
1369MeV 
(154MeV) 4 
1500MeV 



6.0g§ 
130MeV|6] 
782MeV 
0.21 ±0.01 
-0.040 ± 0.003 
0.038 ± 0.003 
0.0017 ± 0.003 
not measured 
0.66fm 
- 1300MeV|7] 



340MeV[|T8] 
510MeV|8] 



>400MeV 



(146 MeV) 4 
~1300MeV 
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